Numerical Strip-Yield Calculation 

of CTOD 
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CTOD background: 
plastic zone sizes 
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■=> Yield zone loading neutralizes stress singularity 
due to remote loading 

■=> plastic zone size estimated from setting K(a+p)=0 


CTOD background: 
plastic zone sizes, cont’d 
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when effects of the crack tip plastic zone are important 


CTOD background: 
some calculation methods 
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Where behind the crack tip to measure CTOD? 

■=> 1 st node, 2 nd node, 45° intercept, or some prescribed distance 
Re-meshing burden for analyses of multiple loads or cracks 


Using Boundary Elements to calculate 
crack face displacements: theory 
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Resulting equation contains hypersingular kernel 

■=> Requires special interpretation; challenging to evaluate 
numerically 


Using Boundary Elements to calculate 
crack face displacements: theory (cont’d) 
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This is the basis of NASGRO’s BE component 


Using Boundary Elements to calculate 
crack face displacements: theory (cont’d) 
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NASBEM is NASGRO’s 
Boundary Element Analysis module 
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NASA and Southwest Research Institute® sign Space Act Agreement 
for joint NASGRO development 

NASGRO Industrial Consortium formed by SwRI; members include 
government agencies and industrial representatives 
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Using NASBEM to calculate CTOD: 

procedure 
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Using NASBEM to calculate CTOD: 
procedure (cont’d) 
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displacement at tip of physical crack a 



Using NASBEM to calculate CTOD: 

results 
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Cracks from holes in infinite sheets 
Periodic cracks in infinite sheet 
3-hole tension specimen 


CTOD verification case: 
reproducing Dugdale’s result 
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CTOD verification case: 
reproducing Dugdale’s result, cont’d 
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CTOD verification case: 
edge crack 
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Verification case: C(T) specimen with W=3, a=1 

• Plastic zone size and CTOD were calculated (p is shown here) 

• Excellent agreement with collocation results by Newman & Mall, 
and also Terada (both 1983) 


CTOD verification case: 

1 crack from a hole, infinite plate 
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CTOD verification case: 

2 cracks from a hole, infinite plate 
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Verification case: plastic zone size studied for various values of a/R 

• Difference between NASBEM and Rich < 2.5% 

• Larger p in small cracks due to higher stress concentration at hole 

• Solution approaches Dugdale solution for large a/R 


CTOD verification case: 
periodic cracks in an infinite sheet 
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• 7 cracks seems a good approximation - idea taken from literature 
on modelling large arrays of fuselage fasteners 

BE compared to Tada (Westergaard stress function, 1974) 


CTOD verification case: 
periodic cracks in an infinite 
sheet, cont’d 
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Relative net section stress, Opet^yield = c/o yield^' a ^) 



CTOD verification case: 
periodic cracks in an infinite sheet, 

cont’d 
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c^b [actual crack + plastic zone / spacing] 



CTOD calculations: 
center crack in finite-width sheet 
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CTOA calculations 
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M(T) specimen: Al 7075-T6, Al 2024-T81 
3-hole tension specimen: Al 7075 


CTOA calculations: 
center crack in finite-width sheet 
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less so for 2024 



CTOA calculations: 
3-hole tension specimen 
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CTOA calculations: 
3-hole tension specimen 
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Summary 
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■=> complex cases with several cracks and holes 
Great accuracy at minimal computational cost 
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Turn capability into production-level tool 

Implement automation of CTOD calculations 

■=> Manual meshing and convergence to K=0 for multiple cracks is 
tedious! 


